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ABSTRACT 



Aims. Scintillation noise is a major limitation of ground base photometric precision. 

Methods. An extensive dataset of stellar scintillation collected at 1 1 astronomical sites world-wide with MASS instruments was used 
to estimate the scintillation noise of large telescopes in the case of fast photometry and traditional long-exposure regime. 
Results. Statistical distributions of the corresponding parameters are given. The scintillation noise is mostly determined by turbu- 
lence and wind in the upper atmosphere and comparable at all sites, with slightly smaller values at Mauna Kea and largest noise at 
Tolonchar in Chile. We show that the classical Young's formula under-estimates the scintillation noise.The temporal variations of 
the scintillation noise are also similar at all sites, showing short-term variability at time scales of 1 - 2 hours and slower variations, 
including marked seasonal trends (stronger scintillation and less clear sky during local winter). Some correlation was found between 
nearby observatories. 
Conclusions. 

Key words. Techniques: photometric - Site testing - Atmospheric effects 



1. Introduction 

One of the main characteristics of astronomical objects is their 
brightness in different spectral bands. The standard precision of 
ground-based photometry is adequate in most cases, but a num- 
ber of astronomical problems require an eve n greater precision 
dHeaslev et aljl996t Everett & Howelll2001l) . One of the funda- 
mental factors limiting the precision of ground-based photom- 
etry is the stellar scintillation occurring in the atmosphere as a 
result of its turbulent nature. 

The fluctuations of the refractive index cause phase distor- 
tion in a plane light wave passing through the atmosphere to 
an entrance aperture of the telescope. As the wave propagates, 
the phase distortions lead to a redistribution of the amplitudes 
between different parts of the wavefront. Averaging within the 
aperture reduces the fluctuations caused by this mechanism, but 
does not eliminate them completely. In photometric practice, this 
effect is considered as an additional source of error. The scintil- 
lation noise expressed in stellar magnitudes does not depend on 
the object's brightness, therefore it cannot be reduced by observ- 
ing brighter stars. 

Scintillation noise in high-precision and fast photometry 
has been studied f or quite a long time (lYoungl 1 19671 119691; 
iDravins et alJI 1997b . The lar ge number of methods proposed for 
reducing scintillation noise dHeaslev et al.l IT996I : IDravins et alJ 
Il998t iGilliland etaTJl 19931; lOsborn et alj|2010h shows that there 
is no perfect solution and that further work is needed. Moreover, 
this same problem is often presented as an argument to do pre- 
cise photometry from space. 

Stellar scintillation is of interest not only to photometry. This 
phenomenon is a powerful tool for remote sensing of optical tur- 
bulence (OT) in the atmosphere. The mechanisms of scintillation 



generation and its characteristics are well studied theoretically 
and experimentally because they are closely related to the most 
important characteristic of the OT above astronomical sites, the 
seeing. 



Unlike the seeing, scintillation received relatively little atten- 
tion in the astro-climatic work; the error budget o f high-precis ion 
photometry is still evaluated using the data from dYounel 19671) or 
other disparate estimates. Only recently such studies have been 
conducted in the general context of the characterization of op- 
tical turbulence ab ove different astron o mical observa tories and 
prospective sites dKenvonetalJl2006l lKornilovll20TTh . These 
data are useful for comparing sites in the context of photometry. 
Moreover, the scintillation noise of long-exposure photometric 
measurements depends not only on the intensity of high-altitude 
OT, but on the wind speed at altitudes above the tropopause 
dKornilovfcOllt) . which is important for the global dynamics of 
the atmosphere. 



This paper presents scintillation noise measurements with 
the MASS instrument at observatories situated in different geo- 
graphical areas. Sect. l2.1l recalls the theoretical description of the 
scintillation in three basic measurement regimes. In Sect. l2.2l the 
method used to estimate relevant parameters from the raw data is 
described. The next Section describes the sites studied here, the 
original data, and the procedure for calculating the scintillation 
noise. Results and comparative analysis are presented in Sect. [4] 
Sect. [5] describes the temporal variability of scintillation noise, 
and the final section is a discussion and comparison with other 
available data. 
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2. Basic relations 

2.1. Theoretical background 

Scintillation is characterized by the variance of relative fluctua- 
tions of the flux / passing through a receiving aperture, the so- 
called scintillation index s 2 : 

s 2 = <(/-</» 2 >/</> 2 . (1) 

In the approximation of weak perturbations, the total scintil- 
lation index is the sum of scintillation produced by independent 
turbulent layers: 

/-»oo 

* 2 = C 2 n (z)W(z)dz, (2) 
Jo 

where W(z) is a weighting function which depends on the size 
and shape of the receiving aperture and does not depend on the 
distribution of the structural refractive-index coefficient C 2 (z). 
The weighting function has a simple physical meaning, it equals 
the scintillation index generated by a layer of unit intensity C 2 dz 
located at a distance z. 

In the general case of non-zero exposure, the weighting func- 
tion also depends on the exposure time and on the wind speed at 
altitude h. The altitude and distance to the layer are trivially re- 
lated as h — z cos y, where y is the zenith angle of the observed 
object. Theoreti cal descr i ptions of the scintillatio n can be found, 
for ex ample, in I Young! (1 19671) : iRoddierl (Il98ll) : iDravins et al.l 
(fl997T) . 

Using the Taylor's frozen-flow hypothesis dTavlorl [l938l) . 
the wavefront evolution during exposure time t for wind speed 
w(z) is reduced to a simple translation by w(z)t. Invoking the 
Kolmogorov spectrum of refractive index per turbations, the fi- 
nal expression for the weighting function is (lTokovininll200l 
lKornilovl201ll) : 

Jr»oo 
&fr m S{z,f)A{f)A s (wr,f). (3) 
o 

Here, the integration is performed over the modulus of the 
spatial frequency /, assuming that the functions in the integrand 
are either axisymmetric or are already averaged over the polar 
angle. The function A(f) is the aperture filter, which takes into 
account spatial averaging by the aperture. A s (w, T,f) is the spec- 
tral filter of the wind translation which describes the temporal 
averaging, and S(z,f) is the Fresnel spectral filter which de- 
scribes the generation of amplitude distortions in the propaga- 
tion of the wavefront. 

For a circular aperture, the aperture filter is A(f) = 
[2Ji{nDf)lnDf] 2 . The Fresnel filter S(z,f) depends on the 
wavelength A of the detected radi ation and can be calculated for 
any spectral energy distribution dTokovininl2003l) . In the case of 
monochromatic radiation S (z, f) = [sin(7T/iz/ 2 )//l] 2 . The charac- 
teristic spatial scale of the filter is the Fresnel radius r? = (Az) 1 ^ 2 . 

The wind-translation filter is not axisymmetric, but after av- 
eraging over the polar angle it is represented as A s (w,r,f) = 
Tiifrw), where the function T\{$) can be expressed in terms 
of special functions dTokovin in 2002; Korn ilovll201 it) . The im- 
portant features of this function are its asymptotes for small £: 
7"i(£) ~ 1 - 7r 2 <f 2 /6 and for large £ Ti (£) * 

It is difficult to analyse the expression (0 in general form be- 
cause it depends on many parameters. However, for certain rela- 
tionships between the parameters the expression is considerably 
simplified. First of all, it is the definition of the small D <K rp 
and large D » rp aperture regimes. In the first case we can set 



A(f) = 1, in the second case we sim plify the Fre snel filter by 
replacing the sine with its argument (Rodd ierlll981l) . 

We can also identify two limiting cases by the value of 
the wind translation. At very short (zero) exposures wt <k 
min(Z), rp), the filter A s (wt, f) = 1 . In the opposite case of long 
exposure wt » max(D, rp), the A s (wt, f) — 1 /jrwrf. 

These limiting situations are generally accepted, although in 
photometric practice, the small aperture approximation is not as 
interesting. In this paper, we will consider the evaluation of scin- 
tillation noise on a large telescope in the regimes of short (SE) 
and long (LE) exposures. 

The scintillation index in a large telescope for sh ort expo- 
sures is expressed by the well-known formula (see, e.g jRoddierl 
EMU): 

s\ = 17.34 ZT 7/3 f C 2 (z)z 2 dz, (4) 
Ja 

since in this case the weighting function is W(z) = 
17.34 D~ 7 t 3 z 2 . The weighting function does not depend on the 
wavelength owing to the achromatism of the scintillation in a 
large telescope. This case is typical for fast photometry. The de- 
pendence on the zenith distance y can be obtained from trivial 
geometric considerations: s 2 , oc M\, where M z = secy. This 
formula ignores the central obscuration of astro nomical tele- 
scopes, important i n the short-exposure regime dYoung|[l967l : 
IDravins et al.lll998l) . 

Nevertheless, in conventional astronomical photometry typi- 
cal exposures of tens of seconds and longer are used. In this case, 
the condition wt » D for long exposure i s satisfied for ty pical 
wind speeds in the upper atmosphere, and dKornilovl2011l) 

st = 10.66 D-^t- 1 [^$fdz. (5) 

Naturally, the wavefront is translated by the wind's compo- 
nent perpendicular to the line of sight. Given that, on average, 
the wind vector is directed horizontally, this effect also depends 
on the air mass M z a nd on the azi muth of the wind with respect 
to the observed star dYoungf T969). The dependence of the scin- 
tillation on the air mass therefore varies from s 2 oc M\ for the 
case of transverse wind up to s 2 oc M\ in the case of longitudi- 
nal wind directed along the object azimuth. In the long-exposure 
regi me the s cintillation is little affected by the central obscura- 
tion ( Yound[l967l) . which has to be ta ken into accou nt only in 
the extreme case of a ring-like aperture dKornilovl2012l) . 

Formulae © and (O show that for a particular telescope and 
exposure time the scintillation index can be easily calculated if 
the integrals i n these formulae ar e known. Using the notation 
introduced in K envon et al.1 d2006): 

S\ = 17.34 f C 2 „(z)z 2 dz = 17.34 M 2 , (6) 
Ja 

for the SE regime, and 

" dz= 10.66 J/ 2 , (7) 

A W (Z) 

for LE regime. 

These scintillation noise parameters S 2 and S 2 character- 
ize the power of scintillation noise under given atmospheric 
conditions and represent the scintillation index in a 1 m tele- 
scope (with 1 s exposure for the LE regime). The quantities M2 
and J/2 are known as a tmospheric moments of second degree 
(Toko vinin etaT]l2003bl) . 
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2.2. Measurement of scintillation noise parameters with 
MASS 

Multi- Aperture Scintillation Sensor (MAS S) is an instrument for 
measu ring OT profile from scintillation dKornilov et al.l [20031 
12007b . Photons from a single bright star collected by a small 
telescope are detected with four photo-multipliers which sam- 
ple intensity in the pupil with four concentric annular apertures 
A, B, C, and D. The smallest circular aperture A at the centre 
has a typical diameter of 2 cm, the largest annular aperture D 
has an outer diameter of ~ 10 cm. These diameters are deter- 
mined by the optical magnification factor k, adjusted and mea- 
sured for each instrument individually. Photon counts with a 1- 
ms exposure are processed statistically to derive 4 normal and 
6 differential scintillation indices which are related to the OT 
through weighting functions (Eq. Apart from the aperture 
geometry, these functions depend on the spectral energy distri- 
bution, i.e. on the spectral type of the star and on the instru- 
ment response. The functions Wih) are computed by numeri- 
cal integration of the basic expression (0. In the calculations, it 
is assumed that the measurement occurs with "zero" exposure, 
neglec t ing the wind translation d uring 1 ms (fTokovini n et al.l 
l2003btlKornilov & Kornilovll201 lb . 

The measured scintillation indices are used not only for the 
restoration of the vertical profile of OT, but also for the eval- 
uation of several atmospheric moments required for comput- 
ing such integral characteristics of the atmosphere as the free 
atmosphere seei ng, isoplanatic angle, e f fective altitude of th e 
turbulence, etc. dTokovinin et al.l l2003bt iKornilov et al.ll2007h . 
Among these moments, the second moment Mi is computed, 
so the evaluation of the parameter S ? is a trivial task. 

The method of finding the parameter S ? has been proposed 
in IKornilov! (1201 lb . It relies on the fact that the scintillation in- 
dex s 2 ^ for long exposure can be calculated from the mean fluxes 
measured with 1 s exposure which are stored in the MASS out- 
put files * . stm. It should be mentioned that for the small MASS 
apertures the LE regime begins at r <K Is. 

In this LE regime, the scintillation index (we call it LE-index 
to distinguish from the "fast" indices with 1 ms exposure) can be 
expressed as 

J A W (Z)T 

where U'iz) are the MASS weighting functions for the LE 
regime. They are distinguished from the conventional weight- 
ing functions by the additional factor {nf)~ l under the integral 
over the modulus of the spatial frequency. Typical be haviour of 
U'iz) f or the whole set of MASS indices is presented in lKornilovl 
(12011b . so we only note that the asymptotic dependencies of 
these functions are U'(z) °c z 4 ^ 3 for infinitely small apertures, 
and U'{z) °c z 2 for large apertures. 

We estimate S ? from the weighted sum (linear combination) 
of LE scintillation indices measured with MASS. The coeffi- 
cients of this sum are adjusted to make the weighted sum of 
U'iz) as close as possible to 10.66 z 2 . Then, as can be seen by 
comparing Eq.|7]and Eq.[8] the weighted sum will approximate 
Sy Details of the procedure are given in appendix lAl 

Fig-fflplots the relative error of such approximations for five 
typical MASS devices. It shows that for propagation distances 
between 4 and 32 km, the errors do not exceed +5%. The contri- 
bution of closest layers (below 1 km) is underestimated, but it is 
negligible due to the factor z 2 . The impact of the distant (beyond 
32 km) OT is underestimated by * 10% as well. Here we ex- 
clude the measurements made far from the zenith (M z > 1.3), so 



1 .2F r 



1.1- 




Fig. 1. The ratio Riz) of the combined LE weighting functions to 
10.66 z 2 . Horizontal lines mark the ±5% and +2.5% corridors. 
The thin solid line depicts this ratio for the D aperture. 



this under-estimation is relevant for altitudes above 25 km where 
the OT has already low intensity. 

This figure also shows the curve for the aperture D of the de- 
vice with k = 16.3. It is evident that if we only use this aperture, 
the result is greatly underestimated. The approximation is there- 
fore better for devices with larger apertures (larger magnification 
coefficient) and a higher spectral response in the blue. 

The influence of the magnification coefficient on the results 
was verified a posteriori for the real data obtained with the de- 
vice MD05 TMT (Armazones). When changing k from 15.9 
to 15.4 (by * 3%), the median value of the parameter S 2 de- 
creased by 4%, so that the requirement for the accuracy of the 
k is not stronger than in the typical processing of MASS data 
(IKornilov et al.ll2007h . The substitution of another spectral re- 
action curve with an effective wavelength smaller by « 14 nm 
changed the median by less than 0.2%. The dependence of the 
approximation on the star colour is negligible. 

3. Data processing 

3.1. Data sets and site testing campaigns 

In the past decade, a large amount of scintillation data from 
MASS instruments has been accumulated. The measurements 
were performed by various projects in different geographical lo- 
cations. Some of these campaigns focused on the site selection 
for future large telescopes. Other studies were methodological 
in nature, or supported already operating telescopes. The list of 
astronomical sites studied in this paper is listed in the Table [T] 

Th ey include all sites studied by the TMT site testing pro- 
gram dSchock et al.ll2009T) . representing two sites in the north- 
ern hemisphere: Mauna Kea (Hawaii) and San Pedro Martir 
(Mexico) and three mountains in Northern Chile: Cerro Tolar, 
Armazones and Tolonchar. The data from this campaign uses a 
similar setup on all sites using MASS/DIMM units attached to 
custom made telescopes capable of operating robotically within 
a large range of wind conditions. 

Measurements on Cerro Armazones in Chile were first car- 
ried out by TMT. After this site was chosen for the E-ELT, moni- 
toring was taken over in 2010 by ESO. Shorter samples at nearby 
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Table 1. Astronomical sites and main characteristics of the site testing campaigns. Observatory coordinates, altitude above sea level 
H and total measurement time T are listed. 
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Fig. 2. Time distribution of the measurements. Indigo colored 
stripes are TMT campaigns, brown stripes are ESO campaigns, 
maroon stripes are CTIO and SAI campaigns. 



sites Ventarrones and La Chira, st udied in the frame o f the E-ELT 
site characterization campaign dVernin et al.l I201 ll) . have been 
added for the purpose temporal correlation analysis. 

The robotic MASS-DIMM site monitor at Cerro Pachon is 
a shared facility between Gemini-S and SOAR telescopes, op- 
erational since 2004. Ea rly results of turbulence monito ring at 
Pachon were reported in lTokovinin & Travouillonl d2006i) . Here 
we use the much larger data set accumulated during 8 years. 

The longest series of observations were obtained at the 
Paranal and Pachon observatories. To control the stability of 
the results and for comparison with the two campaigns at 
Armazones, the Paranal data are also divided in two parts: 2004 
- 2009 (Paranal A) and 2010 - 201 1 (Paranal B). Since the mea- 
surements are carried out at Paranal in continuous mode, both 
samples are quite large. The time coverage of all campaigns is 
shown in Fig. [2] The overall volume of data amounts to more 
than 40000 hours. 

The only site in the eastern hemisphere is represented by 
data from the Mt. S hatdzhatmaz in the Northern Caucasus 
dKornilov et al.l I20T0I) where Sternberg Astronomical Institute 
(SAI) observatory is building its 2.5-m telescope. Note that in 
operating observatories, the measurements with MASS device 
are continuing to provide operational information about OT in 
the atmosphere. 



3.2. Computation of LE scintillation indices 

For each of four MASS apertures, the LE scintillation index 
is calculated as mean square of the difference between adjacent 
1 s mean fluxes: 

s ^27^F§ (F ''- F '' +l)2 - - 001/F ' (9) 

where F is the average of the series of N photon counts {F/} with 
1 s exposure. In the MASS output files, these counts are scaled 
to 1 ms exposure, which causes the factor of 10~ 3 in the last term 
describing the contribution of the photon noise. Use of the differ- 
ence effectively suppresses the contribution of the flux variations 
on time intervals greater than 1 s, unrelated to the scintillation. 
It is easy to show that the scintillation is uncorrelated at long 
exposures; it follows from the relation s 2 (t) oc 1/r in the LE 
regime. 

Each estimate of the uses all 1-s flux values during the 
MASS accumulation time. As a rule, this time is 1 minute, hence 
/V = 60. For an unbiased estimate of the variance of the differ- 
ence, a factor of 1/(N - 1) is required; although the number of 
differences is also N-l, the mathematical expectation is known 
and equal to zero. 

For fluxes registered with 1 s exposure, the photon counting 
statistics can be considered as exactly Poissonian. Any differ- 
ence of the PMTs from an ideal detector does not produce errors 
exceeding 1%, even in the case of extremely small indices. Non- 
linearity of the photon counts is corrected with a fixed dead time 
of 20 ns, as typica l for the MASS d etectors. 

In contrast to iKornilovl d201 lb . the reduction to the zenith 
is performed after evaluation of the parameter § \ from the LE 
indices measured at a certain air mass M z , by Eq. (1A.U . Since 
we do not know the direction of high-altitude winds, two esti- 
mates of the scintillation parameter at zenith are calculated: the 
minimum S ? = S ? M~ 4 for the longitudinal wind, and the maxi- 
mum S? = S\ M~ 3 for the transverse wind (see Sect. 12. U . Data 
obtained at air masses M z > 1 .3 are not used to avoid the uncer- 
tainty in the reduction to the zenith caused by the unknown wind 
direction. In addition, as noted in Sect. 12.21 with a large air mass 
our method underestimates the contribution of high-altitude tur- 
bulence to the parameters S \ and S\. 

Filtering used to remove invalid data and the estimates of 
the resulting errors are considered in Appendix A. The filtering 
procedure eliminates a small, though non-negligible, fraction of 
the data, on average about 10%. In some campaigns, a signifi- 
cant part of measurements was carried out far enough from the 
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Fig. 3. Cumulative distributions of the S \ parameter. The curves 
are labelled from left to right. Four sites are not plotted. 

zenith; in these cases, the proportion of the rejected measure- 
ments exceeds 20%. 

Although we took great care in filtering out the bad data, the 
final results remain almost the same when all the data are used. 
This is so because the data were collected by automatic monitors 
for a long time and therefore contain only a small number of 
faulty measurements. 

4. Results 

4.1. The scintillation parameter S\ 

As noted in Sect. l2.2l the calculation of the parameter S ? is made 
by rescaling the atmospheric moment Aii obtained in conven- 
tional MASS data processing and already reduced to the zenith 
(see Eq. [5). Filtering of the original data consists only in the 
choice of the points corresponding to the correct estimates of 
the parameter S \. For some measurements, for various reasons, 
the AI2 have not been obtained while the corresponding param- 
eter S 3 has been successfully evaluated. These situations are re- 
flected by the fact that in Table [2] the total time T is somewhat 
smaller than in Table [3] 

Note that the algorithm of calculating Mi implemented in 
the early versions of the MASS software had some instabil- 
ity, which caused artifacts in the distribution and a system- 
atic overestimation of the second moment. Therefore, in al- 
most all cases we used the results of data reprocessing carried 
out later with the latest ve rsion of the software atmos-2.93 
dKornilov & Kornilovll201 it) . For measurements on Tololo and 
Pachon we used the Mi calculated from the restored vertical 
profiles of the OT. 

The cumulative distributions of the index S \ (in m 7 ^ 3 ) are 
shown in Fig. [3] and their characteristic points are given in 
Table [2] Of all the distributions, the Mauna Kea site stands out, 
as expected in view of its high altitude and low turbulence in 
the upper atmosphere. The rest of the curves differ only in de- 
tails, especially if we take into account the sampling effects for 
those campaigns with non-uniform seasonal coverage (we dis- 
cuss this in the following Section) and some differences in the 
MASS devices. The fact that the Tololo and Pachon sites show 
a stronger scintillation noise is likely a result of evaluating S\ 



Table 2. Characteristic points of the distributions of the S \ scin- 
tillation parameter. The units are 10~ 4 m 7 ^ 3 . The estimates of rel- 
ative error E and the total accumulated time T are also listed. 
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2.83 


0.06 


2400 


Tolar 


1.29 


1.86 


2.86 


0.09 


1000 


Tololo 


1.48 


2.26 


3.64 


0.10 


1800 


Tolonchar 


1.34 


2.14 


3.64 


0.09 


940 


Ventarrones 


1.44 


2.26 


3.62 


0.09 


1430 



at these sites from the OT profiles. Nevertheless, this systematic 
difference does not exceed 10%. 

The comparison of the Paranal A and Paranal B sub-samples 
shows that the period 2010 - 201 1 does not differ in scintillation 
power from the previous years. The same should be expected for 
the measurements at the Armazones summit due to its geograph- 
ical proximity to the Paranal (30 km). However, the parameter 
S 2 for Armazones B is higher by 15% than for Armazones- A. 
Apparently, this discrepancy is caused by a significant difference 
in the geometry of the entrance apertures of the instruments used 
in both cases. A similar systematic difference is observed for the 
parameter of S3, see Table|3] 

The differential distributions are clearly asymmetric. For all 
the studied sites, the modes of the distributions coincide with 
the 23 - 26% quantiles, so one can use the first quartile from 
Table [2] as the most probable value of S\. The typical skewness 
71 amounts to 2 - 3 and the excess kurtosis is 72 ~ 10. Here we 
do not give more precise values, because adequate calculation of 
the higher moments of a statistical distribution requires careful 
removal of outliers and long-term trends. 

The distributions of In S \ are close to the normal distribu- 
tions with slightly different means and widths. It is also notice- 
able that the asymmetry of the distributions are somewhat dif- 
ferent for different observatories. Typical skewness and excess 
kurtosis of In S ? are small, j\ ~ 0.2 - 0.4 and 72 ~ - 0.2. 

4.2. The 5 3 scintillation parameter 

The relative errors of S 3 estimates are quite large, ~ 0.3 - 0.4, 
widening somewhat the distributions. The significance of this ef- 
fect was estimated by comparing the distributions of the initial 
results obtained with 60 s and with 240 s accumulation time (the 
latter are averages of four consecutive 60-s measures). As ex- 
pected, the distribution of 240-s estimates is somewhat narrower 
than that of the 60-s estimates, but its median is larger by about 
5% because of the asymmetry. 

In Fig. [4] the cumulative distributions of the 240 s estimates 
of S 3 are shown for all sites. For clarity, only the curves cor- 
responding to the assumption of the longitudinal wind (i.e. the 
lower limit of S f) are plotted. In Table[3]the characteristic points 
of the distributions are listed for both wind directions; they can 
be considered as lower and upper limits of the S ? parameter. 
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Table 3. Characteristic points of the distributions of the S ? pa- 
rameter averaged over 4 minutes. The units are 10~ 5 m 4 '' 3 s. The 
estimates of relative error E and the total accumulated time T 
are also listed. 





Long 


itudinal wind 


Transversal wind 


E 


T, h 


5; quartiles 


25% 


50% 


75% 


25% 


50% 


75% 






Armazones A 


0.64 


0.97 


1.48 


0.71 


1.09 


1.66 


0.17 


3820 


Armazones B 


0.73 


1.14 


1.83 


0.82 


1.28 


2.07 


0.17 


1880 


La Chira 


0.82 


1.21 


1.81 


0.92 


1.35 


2.06 


0.16 


530 


Mauna Kea 


0.54 


0.80 


1.23 


0.59 


0.88 


1.34 


0.18 


2270 


Pachon 


0.65 


0.98 


1.49 


0.76 


1.14 


1.74 


0.17 


6800 


Paranal A 


0.70 


1.04 


1.55 


0.80 


1.20 


1.78 


0.16 


5650 


Paranal B 


0.73 


1.08 


1.63 


0.82 


1.22 


1.83 


0.16 


3050 


S. Pedro Martir 


0.76 


1.19 


1.97 


0.82 


1.30 


2.15 


0.16 


1470 


Shatdzhatmaz 


0.60 


0.98 


1.60 


0.66 


1.07 


1.76 


0.16 


2760 


Tolar 


0.78 


1.12 


1.59 


0.89 


1.28 


1.85 


0.16 


1005 


Tololo 


0.65 


0.97 


1.45 


0.76 


1.13 


1.70 


0.17 


1890 


Tolonchar 


0.89 


1.40 


2.25 


1.00 


1.61 


2.59 


0.17 


1000 


Ventarrones 


0.70 


1.07 


1.62 


0.78 


1.20 


1.81 


0.17 


1440 
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Fig. 4. Cumulative distributions of the Si parameter at seven 
sites. The curves for different sites are labelled from left to right. 



Both the curves in Fig. [4] and the data in Table [3] show that 
the distributions of the S \ at various observatories differ even 
less than the distributions of S\. The curve for Mauna Kea shows 
again a lower scintillation. Other observatories, except the SPM, 
have approximately the same characteristics. The scintillation 
parameters measured at SPM are greater than the typical values 
by about 20%. However, we emphasize again that these results 
reflect each particular set of data without taking into account the 
seasonal distribution of the observations. It will be demonstrated 
in Sec. l5.2l that the 5? parameter is subject to seasonal variations. 

The differential distributions of the S\ parameter are very 
similar to the differential distributions for the S\; they have 
roughly the same asymmetry and a somewhat larger excess. Just 
as in the case of S\, the distributions of In S \ are close to nor- 
mal; typical values of their skewness and excess kurtosis are 
n ~ 0.2 - 0.4 and y 2 ~ - 0.5. 



5. Temporal variations of the S 2 and 5 3 parameters 

5.1. Nightly variations 

The study of temporal variations of the astroclimatic parameters 
on time scales of ~ 1 hour (in course of a night) is important 
in terms of s hort-term predictions for operationa l scheduling of 
observations (Racine 1996; Skidmore et al .I2009h . Of course, the 
variability of the parameters S 2 and S 3 is not as critical as the 
variability of the seeing. However, for a long time series of pho- 
tometric measurements, control of the scintillation noise is also 
highly desirable. 

To quantify the nature of the short-term variability, differ- 
ent approaches are used. These are the calculati on of the auto- 
correl ation function (ACF) in the usual sense dTokovinin et al.l 
l2003al) and the average relative or absolute differences of a pa- 
rameter for a certain time d elay. The fractional difference (FD) 
proposed in Racine] (1 19961) is not very amenable to a rigor- 
ous mathematical analysis (being a ratio, it has a non-trivial 
distribution function) . The absolute difference (A D) used in 
Skid more et alj d2009l) ; iGarcfa-Lorenzo et al.l (1201 Ol) is restrict- 
edly suitable for non-stationary processes where the variance de- 
pends strongly on the average. 

To achieve the objective outlined by Racine1 (ll996l) in his in- 
troduction of the FD, one could analyze logarithms of the param- 
eters instead o f the parameters themsel ves. This method is rec- 
ommended by iJenkins & Watts! (1 19691) for non-stationary pro- 
cesses of multiplicative nature, where the dispersion is nearly 
a linear function of the mean. We have already used such a 
transformation when evaluating random errors of S \ and S\. In 
Sect. [4] we pointed out that In S\ and In S ? are distributed almost 
normally. 

Thus, to describe the temporal behaviour of the scintillation 
parameters we use the structure functions of their logarithms. 
For the S \ parameter: 



D 2 (T) = {(\nSl(t)-\nSl(t + T)f) t , 



(10) 



where the averaging is performed over all available pairs for the 
given delay (lag) T. Similarly, we define the function Dt,(T) for 
the parameter Si. In order not to complicate the analysis with 
normalization, we consider non-normalized functions and define 
the characteristic correlation time of the scintillation power us- 
ing an absolute criterion. 

The structure functions for all sites are presented in Fig. 
and Fig. [6] where the contribution of uncorrected noise was sub- 
tracted. For this calculation we used 1 minute measurements of 
both S j and S \. Note that all structure functions have a simi- 
lar behaviour for delays less than 2 hours. At longer delays, the 
time coverage of the individual data sets influences the computed 
structure functions and their behaviour becomes irregular. The 
initial part of the functions 'D^iJ) and T>^(T) is well approxi- 
mated by the formula 



D(T) 



■exp(-y/f/Tc + S)]. 



(11) 



The main features of this dependence are infinite derivative 
at zero and saturation at the level of <x 2 when T — > 00. The point 
T = t c corresponds to 0.632 £)(oo) = 0.632 <x 2 . The time con- 
stant t c therefore characterizes relative changes in the power of 
the scintillation. The structure function (TiTT i corresponds to the 
temporal spectral density proportional to oc 1/cd (flicker noise) 
which saturates at low frequencies o> <k \jnr c . 

The approximation of the measured structure functions by 
the Eq. (fTTT i on the initial segment T < 2 hours shows that the 
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Fig. 5. Structure function £>2(T) of logarithm of the scintillation 
parameter S\ at seven sites. Thin solid lines depict the approxi- 
mations by Eq. [TTI 
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Fig. 6. Structure function Di(T) of logarithm of the scintillation 
parameter S 2 at seven sites. Thin solid lines indicate the approx- 
imations by Eq.fTTl 

parameters cr 2 and t c are similar at all sites, with typical val- 
ues t c ~ 1 hour. (Except for the Tolonchar site, where the time 
constant is longer than 10 hours). The values of r c are given in 
Table [4] and [5] The parameter a 2 lies in the range of 0.25 - 0.4. 
The small additional parameter 5 ~ 0.05 compensates for the 
imperfect subtraction of uncorrected noise, it is not considered 
here. 

For practical purpose, another way to evaluate the stability 
of the scintillation noise may be more convenient. To do this, we 
transform the difference of the logarithms in ( TTOb to a logarithm 
of ratio and define the statistics 

R 2 (T) = exp(D 2 (T) 1/2 ), (12) 

and the similar function R?,(T) for the S 2 parameter. The func- 
tions Ri(T) are plotted in Fig. [7] The time lags T1.5 corresponding 
to changes of S 2 and S 2 by 1.5 times are given in Table |4] and 




1.1- 



•1 ni i i i 1 i i i i i i i i i i i 

1/32 1/16 1/8 1/4 1/2 1 2 4 

T, h 

Fig. 7. Ratio of S 2 parameters for a time lag T, Rj(T), at seven 
sites. The ratio of 1.5 times is shown by the horizontal line. 

At almost all sites, the S 2 parameter is changing faster than 
the S\. Apparently, in this case the variations of the wind speed 
are important. The fastest changes of the scintillation power are 
observed at Mauna Kea: 0.2 and 0.3 hours. The atmosphere 
above Mt. Shatdzhatmaz is changing most slowly: 0.7 and 0.8 
hours. However, note that most measurements at this site have 
been performed in the fall which is the most stable season. 

Variation of the scintillation power over longer time scales 
can be estimated by comparing its variability within a night to 
the variability during the measurement season. To do this for 
both scintillation parameters, we calculated the average (In S 2 )j 
and the corresponding variances cr 2 within each night j. Then, 

the variance of the average nightly quantities cr 2 was evalu- 
ated. The intra-night variability is quantified by its mean over 
the whole campaign, {cr 2 .). Nights containing less than 10 mea- 
surements of the S 2 and S 2 were discarded in this analysis. 

The (cr 2 ), corrected for the contribution of random measure- 
ment errors, are presented in Table|4]and|5] On average, the vari- 
ance of the random noise is 0.01 for the S \ parameter, and about 
0.1 for the S 2 (see Table [2] and 0. It is seen that the relative 
variability within a night is almost identical at all observatories. 
As expected, the variability of the In S 2 is slightly higher than 
the variability of the In S\. There is one statistically insignificant 
exception for the measurements carried out at Tololo. The aver- 
age value of the variability within nights for \nS\ is 0.13 ± 0.01, 
for lnS 2 is 0.16 ±0.03. 

The variability of the nightly averages is a combination of 
regular seasonal changes (see the next Section) and the spo- 
radic variability from night to night. We did not isolate the sea- 
sonal behaviour, so the full variances cr 2 , are listed in Table [4] 
and [5] Not surprisingly, the maximum values are observed for 
Shatdzhatmaz and S. Pedro Martir, as the seasonal variability 
at these sites is as large as 0.1 - 0.14 (about 2.5 times larger 
than the intra-night variability). Measurements at Tolonchar also 
show a significant long-term variability. Many parameters de- 
scribing the scintillation noise at this site are substantially dif- 
ferent from the other sites. The median variability from night to 
night at Tolonchar is 0.31 + 0.07 for In S 2 and 0.24 + 0.08 for 
InS 2 (less than for In S 2 , which is atypical). 
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Fig. 8. Seasonal variation of the scintillation parameter S\. Monthly medians are plotted as circles, the quartiles are indicated by 
bars. The grey bars show the share of the measurements (right axis) entering the respective month. 



Table 4. Variation of In S \ during the night and between the 
nights. The time constants r e and T1.5 are given in hours. 



Table 5. Variation of In S 2 during the night and between the 
nights. The time constants t c and t\ .5 are given in hours. 



Site 


N„ 


T 


°i 






T1.5 


Site 


N„ 




< 


«r}> 




n.5 


Armazones A 


928 


0.399 


0.267 


0.130 


1.53 


0.43 


Armazones A 


929 


0.487 


0.206 


0.173 


0.82 


0.26 


Armazones B 


525 


0.459 


0.310 


0.132 






Armazones B 


526 


0.559 


0.244 


0.215 






La Chira 


201 


0.411 


0.271 


0.124 


2.46 


0.36 


La Chira 


202 


0.495 


0.240 


0.191 


0.82 


0.18 


Mauna Kea 


576 


0.498 


0.361 


0.141 


0.58 


0.31 


Mauna Kea 


577 


0.501 


0.232 


0.169 


0.48 


0.21 


Pachon 


1426 


0.469 


0.339 


0.144 


1.79 


0.33 


Pachon 


1770 


0.488 


0.243 


0.149 


2.36 


0.37 


Paranal A 


1365 


0.416 


0.293 


0.123 


1.01 


0.46 


Paranal A 


1371 


0.465 


0.229 


0.156 


1.43 


0.29 


Paranal B 


582 


0.406 


0.281 


0.125 


1.01 


0.46 


Paranal B 


583 


0.457 


0.216 


0.179 


1.43 


0.29 


S. Pedro Martir 


608 


0.519 


0.378 


0.126 


2.08 


0.48 


S. Pedro Martir 


609 


0.570 


0.348 


0.144 


2.00 


0.32 


Shatdzhatmaz 


706 


0.472 


0.417 


0.094 


0.84 


0.80 


Shatdzhatmaz 


712 


0.621 


0.451 


0.106 


0.97 


0.68 


Tolar 


264 


0.355 


0.209 


0.119 


1.04 


0.55 


Tolar 


264 


0.388 


0.149 


0.116 


0.53 


0.33 


Tololo 


625 


0.427 


0.328 


0.121 


1.13 


0.42 


Tololo 


630 


0.460 


0.258 


0.113 


1.00 


0.45 


Tolonchar 


386 


0.575 


0.477 


0.135 


>10 


0.44 


Tolonchar 


409 


0.619 


0.342 


0.176 


>10 


0.27 


Ventarrones 


553 


0.439 


0.299 


0.124 


0.99 


0.47 


Ventarrones 


563 


0.486 


0.234 


0.154 


0.69 


0.27 



5.2. Seasonal variability 

For the all observatories, the long-term evolution of the param- 
eters S 2 and S \ is seasonal. Seasonal variability is evident in 
the raw data, but for a more detailed and quantitative study, we 
calculated the statistical characteristics by months of the year. 



These medians and quartiles (in the form of bars) for the param- 
eter S 2 are shown in Fig. [8] and for S \ in Fig. [9] 

Note that almost all the samples on which these statistics 
were determined are sufficiently large and their volumes vary 
by no more than 2 times from month to month. The exception 
is the Mt. Shatdzhatmaz data, where the vast majority of clear 
weather happens during the fall. The seasonal dependence of S 3 
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for Armazones B (not plotted) is similar to that of Armazones A, 
with slightly larger values. 

In these plots, the noticeable difference between the northern 
and southern hemispheres is seen, especially for 5 1 ?. There is 
also some dependence of the seasonal variability on the latitude 
of the observatory. The lowest seasonal variability is found at 
Mauna Kea (for both S \ and S f). The most pronounced seasonal 
variability of S \ can be seen on the curve for Shatdzhatmaz in the 
summer months; the scintillation noise increased by almost two 
times, although this makes little effect on the overall median (see 
Table [3}. In contrast, the data for the SPM observatory contain 
a large number of measurements from May to August, which 
over-estimates the median for the entire data set. 

Conversely, in the southern hemisphere the median of S ? is 
minimal in the period from May to September, and the median 
S 2 is maximal from July to September, while in the northern 
hemisphere there is a minimum in this season. However, in gen- 
eral the seasonal variability does not exceed the amplitude of dis- 
tributions defined by the quartiles. This means that at any time 
of the year there is a reasonable probability of both good and bad 
conditions for high-precision photometry. 

Seasonal variations in the scintillation noise are caused by 
the redistribution of seasonal winds in the upper atmosphere. 
Moreover, an increase in the wind speed leads to an increase 
in the S 2 parameter but also to a decrease of the S\. Data on the 
wind speed above the tropopause may also be obtained from our 
measurements. 

5.3. Correlation between close observatories 

Strong correlation of scintillation power on time scales of the or- 
der of 0.5 - 1 hours (Sect. 15. U implies that the size of the spatial 
region in which the power of turbulence, wind, and their distri- 
bution over altitude can be considered constant, is ~100 km. In 
such case, a strong correlation of the scintillation power should 
be observed for closely spaced observatories (<50 km). 

Some observatories on our list are geographically located 
very close to each other. The distance between Tololo and 
Pachon is about 10 km. The distance between Paranal and 
Armazones is also no more than 25 km. Additionally, in an area 
around Paranal, the ESO team investigated two sites: La Chira 
(15 km) and Ventarrones (32 km). Two other summits in north- 
ern Chile were studied by the TMT team: Tolar and Tolonchar 
are spaced from Paranal by about 300 km, first to the North and 
the other to the West. 

The cross-correlation of the scintillation parameter S \ (or 
Sf) measured at two observatories, was calculated as follows. 
For the required time delay T, the pairs of measurements for 
which the difference between their acquisition time is less than 
T/2 (30 s) were selected from the two data sets. Then, using 
the classical method, the correlation coefficient p(T) of the first 
sample on the second for the logarithms of the scintillation pa- 
rameters In S 2 was calculated. Finally, we correct the correlation 
coefficient for the contribution of random measurement errors by 
its renormalization with corrected variances. 

Naturally, the correlation reflects not only the short-term 
variations of the scintillation power, but also its variability on 
the time scale of days, including seasonal variability. This can 
be clearly seen in Fig. [TOl which shows the cross-correlation be- 
tween Tololo and Pachon with time lags -5 < T < +5 hours. 
In addition to the correlation peak of ~ 1 hour width caused by 
the intra-night variability, there is a wide pedestal at the level 
of p(T) * 0.5. The cross-correlation between measurements at 




Fig. 10. Cross-correlation between the In S \ measurements at 
Tololo and Pachon. For comparison, the auto-correlation func- 
tion for Tololo is plotted (dashed line). The cross-correlation be- 
tween Tololo and Paranal (crosses) is close to zero. 



Tololo and Paranal, also depicted in the figure, shows no cor- 
relation peak, while its constant level of p(T) w 0.2 reflects 
the common nature of seasonal variability at observatories in 
Chile. In contrast, a typical cross-correlation between Paranal 
and Shatdzhatmaz is « -0.2 (the seasonal variability is anti- 
correlated). 

The width of the correlation peak is comparable to the width 
of the peak of the auto-correlation function (ACF) for each in- 
dividual site. For illustration, the ACF for Tololo is plotted in 
Fig. [10] A similar situation is observed for other close observa- 
tories, see Fig. QT| except for weak cross-correlation peak be- 
tween Paranal and Ventarrones. The maximum cross-correlation 
is larger for the parameter S \ than for S 2 , but not by a significant 
amount. The height of the peak above the extended pedestal is 
approximately 0.2. Generally, when there is a significant corre- 
lation, the ratio of the cross-correlation peak to the pedestal is in 
good agreement with the relative variability during the night and 
the variability from night to night (Tables [4] and [5]). 

The lack of a significant delay in the variations of the scintil- 
lation power indicates that at large spatial scales (10 - 100 km), 
the wind transportation of the OT is not dominant compared to 
the intrinsic evolution of the turbulence. 



6. Discussion and conclusion 

The resulting characteristics of the scintillation noise in the 
regime of short exposures are shown in Table [6] for long ex- 
posures - in Table [7] Instead of the squared parameters, those 
tables list S 2 and S 3. The S 3 was calculated as the average be- 
tween the values obtained under the assumptions of longitudinal 
and transverse wind. For the Armazones site, the average val- 
ues computed with a weight proportional to the number of mea- 
surements for campaigns Armazones A and Armazones B are 
presented. 

Recall that the scintillation noise at the zenith cr s and <x L for 
a telescope with diameter D is calculated by the formulae 

CT S =S 2 D- 1/b (13) 
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Fig. 9. Seasonal variation of the scintillation parameter S 2 . Monthly medians are plotted as circles, the quartiles are indicated by 
bars. The grey bars show the share of the measurements (right axis) entering the respective month. 



for short exposures, and 

ct l =S 3 D- 2I3 t- 112 



(14) 



for a long exposure r. 

These expressions follow directly from the Eqs. and ©. 
To convert the values of the Table[6]and[7]into stellar magnitudes, 
they should be multiplied by the constant 1.086, i.e., the ampli- 
tude of the scintillation noise in magnitudes is <x mag = 1 .086 cr. 



As we pointed in Sect. 12.11 for the SE regime the central 
obscuration ef fect should be c onsidered. A good approximation 
is presented in Kor nilovl (120121) . For extra-large telescopes (D > 
8 m), the effect of the turbulence outer scale is substantial in both 
the SE and LE regimes, and should also be taken into account 
dKornilovl2012l) . 

The results shown in these tables are in good agreement with 
the estimates of scin tillation parameters S 2 and S 3 for Tololo and 
Pachon obtained in iKenvon et al.l (120061) by numerical calcula- 
tion of the moments M2 and if 2 by Eqs. (|6]l and O on the basis 
of the measured OT vertical profiles and the modelled wind pro- 
file. Our method gives fo r the Maidanak observatory a similar 
estimate of S3 = 0.0030 dKornilovll201 lb . Unfortunately, other 
statistically reliable estimates of the scintillation noise are not 
available to our knowledge. 

The results indicate that scintillation noise is mostly de- 
fined by the global turbulence at altitudes of 10-15 km, at the 



tropopause and above. The only site in our list where the lo- 
cal effects are noticeable and make a difference with other sites 
is the Tolonchar. In this case, the proximity to the main ridge of 
the Andes which is quite high and perpendicular to the global air 
circulation could create a quasi-stationary vertical vortex with a 
scale of tens of kilometres. 

The stronger scintillation noise (about 10% excess) at the 
SPM observatory is likely a result of a biased seasonal distri- 
bution of the observations together with a significant variability 
from night to night (see Table [5}. On the other hand, the most 
probable value of S 3 at this observatory exceeds the value at 
Paranal by only 3%. 

6.1. Comparison of the S3 parameter with the Young's 
equation 

Until now, astro nomers (see, e.g., lEverett & Howelll 2001: 
iMann et all 1201 ll) used the Young's formula (lYoung 19671) to 
estimates the contribution of scintillation noise to the accuracy 
of photometric measurements: 



0.0030 D- 2/3 M 3 ' 2 e-WAo-r-i/a 



(15) 



where D is telescope diameter in meters, and /! b s is the observa- 
tory altitude above sea level. Since the original expression con- 
tains a bandwidth rather than exposure, there was a misunder- 
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Fig. 11. Same as in Fig. [10] for Paranal and the other nearby 
sites. For comparison, the auto-correlation function for Paranal 
is plotted in dashed line. Short-term cross-correlations between 
Paranal and Tolar (crosses) and Tolonchar (pluses) are absent. 

Table 6. Comparison of the S 2 distribution at all sites (the units 
are m 7 ^ 6 ). 







52 quartiles 






25% 


50% 


75% 


Armazones 


0.0108 


0.0132 


0.0166 


La Chira 


0.0127 


0.0155 


0.0198 


Mauna Kea 


0.0088 


0.0109 


0.0139 


Pachon 


0.0124 


0.0155 


0.0197 


Paranal 


0.0117 


0.0144 


0.0181 


S. Pedro Martir 


0.0117 


0.0147 


0.0192 


Shatdzhatmaz 


0.0107 


0.0132 


0.0168 


Tolar 


0.0113 


0.0136 


0.0169 


Tololo 


0.0122 


0.0150 


0.0191 


Tolonchar 


0.0116 


0.0146 


0.0191 


Ventarrones 


0.0120 


0.0150 


0.0190 



Table 7. Comparison of the S 3 distribution at all sites (the units 





5 3 Mode 


25% 


Si quartiles 
50% 


75% 


Armazones 


0.00300 


0.00266 


0.00329 


0.00410 


La Chira 


0.00328 


0.00294 


0.00357 


0.00439 


Mauna Kea 


0.00265 


0.00238 


0.00290 


0.00358 


Pachon 


0.00298 


0.00265 


0.00325 


0.00402 


Paranal 


0.00311 


0.00275 


0.00336 


0.00410 


S. Pedro Martir 


0.00310 


0.00280 


0.00352 


0.00453 


Shatdzhatmaz 


0.00280 


0.00251 


0.00320 


0.00409 


Tolar 


0.00329 


0.00290 


0.00346 


0.00414 


Tololo 


0.00303 


0.00265 


0.00324 


0.00396 


Tolonchar 


0.00353 


0.00306 


0.00386 


0.00488 


Ventarrones 


0.00320 


0.00272 


0.00337 


0.00413 



standi ng in the translation of one to another. The (fT~5T > is taken 
from (Gilliland et al.l 19931) . where it was corrected after the in- 
tervention of Young ( 19931) . The d ependence on t he observatory 
altitude was originally proposed by Reiger ( 1963), who assumed 
an exponential dependence C\ from altitude, but the scale height 
of ho = 8 km was established by Young. The numerical coeffi- 



cient was determined from observations mainly with the 0.9 m 
telescope. 

For comparison with our data, we rewrite the previous for- 
mula as S3 = 0.0030 e~ h ° h *l ha . Using this formula and the site al- 
titudes from TableQ] we obtain S3 estimates from 0.0018 for the 
highest summit of Mauna Kea to 0.0023 for the Shatdzhatmaz. 
A comparison with Table [7] shows that the Young's formula 
underestimates the median amplitude of the scintillation noise 
by a factor of 1.5 (scintillation power by two times). The val- 
ues of S 3 inferred from the Young's formula are similar to the 
first quartiles of its actual distribution s derived both here and in 
iKenvon et al] 120061) : iKornilovl (T20T1 . 

The dependence of the scintillation on the observatory alti- 
tude is ambiguous because of many factors affecting it. Our re- 
sults do not show this dependence for the most of studied sites 
located between 2000 and 3000 m a.s.l. because this effect does 
not exceed the accuracy of the method. Only for the highest site, 
Mauna Kea, we can possibly relate lower scintillation noise to 
higher altitude. On the other hand, the Mt. Ventarrones data show 
domination of local effects. 



6.2. Conclusions 

This paper presents the results of the evaluation of scintilla- 
tion noise in observations on telescopes with large diameter 
D » rp or D > 1 m for the optical and near-infrared in the 
Earth's atmosphere. These estimates are obtained by an indi- 
rect method based on the data of the measurement with the 
MASS instruments, without involving measurements on large 
telescopes. However, 1) this method has a reliable theoretical 
basis, 2) identical instruments were involved and measurements 
were obtained using the same technique, 3) data for each stud- 
ied site were obtained over a long time period and have a large 
statistics. 

The scintillation nois e at short exposures i s quantified by the 
parameter S 2 proposed in lKenvon et alj d2006l) . while the param- 
eter S 3 characterizes observations with long exposures. Using 
these parameters, one can calculate the scintillation noise for 
a telescope of any reasonable diameter for any long exposure. 
For example, 8-m telescopes such as Gemini or LSST at Cerro 
Pachon will have scintillation noise of 50 //mag for r = 5 min on 
a typical night. 

The data from seven of the 1 1 sites studied here should be 
of particular interest to astronomers because they describe the 
conditions at the existing observatories Tololo, Pachon, Paranal, 
S. Pedro Martir, Mauna Kea, or the observatories soon to come 
into operation (Armazones, Shatdzhatmaz). Scintillation noise 
measurements at four other summits in the northern Chile are 
more interesting from the point of view of general turbulence 
behaviour in the upper atmosphere. 

In addition to the general statistical characteristics of the 
scintillation noise parameters, their temporal variability was in- 
vestigated, from the short time scales of minutes to seasonal 
variations. Characteristics of the short-term (during the night) 
changes are important for optimizing observational strategy in 
the high-precision photometry. Usually we can assume that the 
scintillation noise power is stable enough for a 0.5 - 1 h, but 
sometimes it may change significantly on timescales of 10 - 
15 min. 

The seasonal variations are significant enough and in some 
observatories they reach a factor of two in power. They are di- 
rectly related to the latitude of the observatory. The minimum 
amplitude is observed for Mauna Kea which is close to the 
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equator and the maximum for the mid-latitude observatory at 
Mt. Shatdzhatmaz. 

The main conclusion from the comparison of scintillation 
noise at different observatories is that there are no major dif- 
ferences. For the purpose of choosing the best site for high- 
precision photometry, all sites are essentially equal. A much 
larger effect can be achieved by choosing the best season for 
observations and by using the real-time information about the 
power of the scintillation in operational planning of the photo- 
metric observations. 

Apart from the high-precision photometry, the parameters 
measured in the paper can be interesting for erro r budget eval- 
uation of high-precision d ifferential astrometry dKenvon et alJ 
2006: ICameron et al.ll2009h . 
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Appendix A: The weighting function approximation 

The best approximation is found in the same way as for other at- 
mospheric moments dTokovinin et al.ll2003bl) . We solve a linear 
system of equations U'd = q, where U' is the weighting function 
matrix of dimension m x n, d is the vector of unknown coeffi- 
cients, q is the desired function 10.66z 2 . The number of nodes 



Table A.l. Coefficients dj of the approximation of the function 
10.66 z 2 by a sum of weighting functions U'(li) for the different 
MASS/DIMM devices. The devices are marked with their as- 
sembly number and project. The magnification coefficient k and 
noise factor F^ are listed as well. 



Site 


Device 


k 




da 


dc 


d D 




Armazones A 


MD05 


15.9 


-0.009 


0.008 


-0.029 


0.107 


1.47 


Armazones B 


MD31 


14.0 


-0.014 


0.020 


-0.054 


0.117 


1.89 


La Chira 


MD21 


16.1 


-0.009 


0.006 


-0.024 


0.104 


1.40 


Mauna Kea 


MD08 


15.8 


-0.010 


0.008 


-0.031 


0.108 


1.48 


Pachon 


MD07 


14.5 


-0.011 


0.015 


-0.047 


0.115 


1.74 


Paranal 


LITE 


16.8 


-0.009 


0.004 


-0.020 


0.103 


1.37 


S. Pedro Martir 


MD11 


15.8 


-0.009 


0.006 


-0.024 


0.101 


1.40 


Shatdzhatmaz 


MD09 


16.3 


-0.009 


0.007 


-0.026 


0.107 


1.42 


Tolar 


MD05 


15.9 


-0.009 


0.008 


-0.029 


0.107 


1.47 


Tololo 


MD02 


15.0 


-0.010 


0.012 


-0.040 


0.112 


1.63 


Tolonchar 


MD02 


15.0 


-0.010 


0.012 


-0.040 


0.112 


1.63 


Ventarrones 


MD21 


16.1 


-0.009 


0.006 


-0.024 


0.104 


1.40 



n = 50 of the distance grid {z,} is substantially greater than the 
maximum number of indices m = 10. As always, a log-uniform 
distance grid is used, with higher density at low altitudes. 

The system is solved by singular value decomposition, dis- 
carding singular values less than 5 • 10~ 4 . The quality of the 
approximation is controlled by the noise amplification factor 
f N = (Hjd 2 ) 1 * 2 / Yijdj- This system is weighted oc z~ l to en- 
sure best fit in the range of 5 - 25 km. Here we use only the 
four normal indices because it became clear that an approxima- 
tion which includes the differential indices is better only at low 
altitudes, the least interesting zone from the standpoint of the 
scintillation in large telescopes. 

Using the linear relationship between the weighting func- 
tions and the scintillation indices, we can finally write 

S\ = Yi d &P 7 = A,B,C,D, (A.l) 

j 

where s? . are the LE indices (with 1 s exposure) measured in the 
MASS apertures A, B, C, and D. 

The weighting functions U'(z) and, therefore, the coefficients 
dj, depend on the aperture dimensions and spectral sensitivity 
of the MASS detectors. Although all MASS/DIMM devices are 
almost identical, they have been used with different feeding op- 
tics and somewhat differ in the spectral response of the detec- 
tors and in the magnification coefficient k for overall optical sys- 
tem "MASS+telescope" that determines the aperture diameters 
in plane of the entrance pupil. Therefore, the coefficients dj were 
computed for each instrument considered in this work individu- 
ally. These values are given in Table I A. ll 

Appendix B: Filtering of the LE scintillation indices 

During the computation of the indices using the Eq. (0 and the 
data stored in the file * . stm, a set of the flux values {F,} is omit- 
ted if 1) there is no information about the sky background or 
the background is too high, 2) the flux ratio between MASS 
apertures B and C does not satisfy the condition Fr < Fq, 3) 
the number of 1-s flux values in a 1-min. accumulation period 
TV < 10, and 4) the flux in aperture D is low, Fu < 100. These 
criteria eliminate obviously wrong or inaccurate data. However, 
they do not exclude all situations with incorrect data and we de- 
veloped additional data-filtering criteria: 
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5 10 15 20 25 30 



g 2 

Fig. B.l. Distribution of the g2 indicator for all data from 
Shatdzhatmaz site. The dashed line represents the complemen- 
tary cumulative distribution. The thin vertical line shows the 
chosen cut-off threshold. 

- g\ is the ratio of the maximum term in eq. (0 to the total sum 
for aperture D. This parameter detects strong outliers. 

- g2 the ratio of the flux range in the aperture D between 
its 95% and 5% quantiles (i.e. ignoring the 2 lowest and 2 
largest values in series of 60 points) to its rms fluctuations 
Vs 2 . This indicator is sensitive to a trend. 

- g3 is calculated similarly g2 using the full flux range. It mon- 
itors isolated overshoots or, in their absence, the general 
trend. 

We determined the filtering thresholds for these indicators from 
their empirical distributions. 

As an example, the differential and complementary cumula- 
tive distributions of g2 for the data obtained at the Shatdzhatmaz 
are shown in Fig. IB. II Similar distributions are found at other 
sites. The inflection point on the right (descending) branch of 
the differential curve is chosen as the threshold for all sites (at 
larger g2 the distributions clearly change their character). In this 
example, the cut-off is set at g2 > 6.7, which discards m 8% of 
the data. Visual inspection confirms that the discarded data are 
affected by either clouds or large tracking errors. 

Using differences for calculating the LE scintillation indices 
does not completely suppresses the low-frequency flux variation, 
increasing the estimated variance. Most of these variations have 
the form of a trend during the accumulation time. We denote the 
relative change in the flux due to the trend by A. It is easy to 
show that a trend adds A 2 /2N 2 to the true value of the LE index 
s 2 ^. If the relative error of the index estimate must not exceed p, 
the following condition must be satisfied: 

A < N(2p s 2 ) 112 . (B.l) 

If p — 0.01 and N = 60, the allowable trend must be 
A < 8.5 {s 2 ) 1 ^ 2 . Since the LE indices in the apertures A, B and 
C enter in the result (the parameter S 2 ) with smaller coefficients, 
and are themselves larger than the index in D, it is sufficient to 
check the condition ( IB. Il l for the aperture D only. In addition, 
monitoring of the absolute value of the trend is also needed, 
since large systematic changes in flux induce such an increase 




1 2 3 4 5 



Fig. B.2. Distributions of the ratios of the s£ indices in the aper- 
tures A, B, C to the index in D measured at Mauna Kea. Vertical 
thin lines indicate the theoretical ratio for the distance z — 2 km 
(top) and 32 km (bottom) for the C/D, B/D and A/D, from left to 
right. 



in the variance that the condition ( IB. II ) again begins to be ful- 
filled. 

We used yet another powerful criterion based on the theoret- 
ical relation between the measured values. The behaviour of the 
weighting functions U'(z) restricts the ratios of the LE indices. 
Clearly, that the ratios U'(z) for the apertures A, B, and C, to 
the weighting function for the aperture D not equal to the ratio 
of the indices themselves, but for any profiles of OT and wind 
the ratios of s 2 ^ must lay between the theoretical minimum and 
maximum. 

The ratio of LE indices Jlc/ s ld' s weakly dependent on the 
distance, which is reflected in the form of the observed distribu- 
tion of the ratio of in the apertures C and D shown in Fig. IB.2l 
The cumulative distribution is very steep and the differential dis- 
tribution is very narrow. The ratio of indices in the apertures A 
and D typically has a broader distribution, which is explained 
by the behaviour of the functions U'(z). This ratio is most sensi- 
tive to the specific characteristics of the MASS device and to the 
minimum altitude of the turbulence which produces noticeable 
scintillation. The minimum ratio of LE indices is reached when 
all the turbulence is at ~ 32 km. 

Fig. lB.2l indicates a good match between the measured ratios 
and their theoretical estimates, so we boldly used the experimen- 
tal distributions to determine the lower and upper limits (approx- 
imately at the level of 0.1%) for rejecting the outliers. Typical 
instrumental causes of such outliers are a wrong sky background 
value or an accidental partial vignetting of the entrance pupil. Of 
course, random errors in the LE indices also widen the observed 
distribution. 



Appendix C: Error estimation of LE indices 

The errors of LE indices s? were estimated in two ways. First of 
all, assuming stationarity over one minute, normal distribution, 
and uncorrelated 1-s fluxes, we estimated the probable relative 
error of the sample variance s = (2/N) 1 ^ 2 as 0.18 for all aper- 
tures. 
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The second method consists in computing the mean squared 
difference between the logarithms of successive estimates: 5, = 
(In sf - In s 2 +1 ) 2 . Assuming a quasi-stationarity, the relative error 
is e — (0.5(c))) 1 / 2 . The probability distribution of 5 is similar to 
the x distribution with one degree of freedom. The medians of 
such estimated errors are thus 0.175, 0.174, 0.172 and 0.172 for 
the apertures A, B, C and D (for the Shatdzhatmaz data). 

Assuming the same relative errors of s? for all apertures, 
s = 0.18, and assuming uncorrected errors in all apertures, 
we calculated the relative error E\ of the parameter S 2 . Because 
noise amplification by the linear combination ( IA.U depends on 
the coefficients which differ between individual MASS instru- 
ments, the errors of S 2 are also different. The median values of 
these errors are 0.28 - 0.36. 

The errors in the different apertures are partially correlated 
(see lKornilov & Kornilovl201 II) . so the estimate E\ of the errors 
was additionally checked by calculating the mean difference of 
the logarithms of the adjacent S 2 values, as for the LE indices. 
These estimates Ej are somewhat smaller than E\, but there is a 
good agreement between them. 

The average errors of the 4-minute S 2 estimates obtained in 
the this way are listed in Table [3] It can be seen that they are 
practically identical for all sites. The same method was used 
to estimate the errors of the parameter S\, which are typically 
0.06 - 0.09 (see Table 0. These errors are smaller because the 
estimates are calculated from ~ 1000 samples. 



14 



